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Abstract Wc present two complementary approaches for determining the ref- 
erence for the covariant Hamiltonian boundary term quasi-local energy and 
test them on spherically symmetric spacetimes. On the one hand, we isomct- 
rically match the 2-surface and extremize the energy. This can be done in two 
ways, which we call programs I (without constraint) and II (with additional 
constraints). On the other hand, we match the orthonormal 4- frames of the 
dynamic and the reference spacetimes. Then, if we further specify the observer 
by requiring the reference displacement to be the timelike Killing vector of the 
reference, the result is the same as program I, and the energy can be posi- 
tive, zero, or even negative. If, instead, we require that the Lie derivatives of 
the two-area along the displacement vector in both the dynamic and reference 
spacetimes to be the same, the result is the same as program II, and it satisfies 
the usual criteria: the energies are non-negative and vanish only for Minkowski 
(or anti-de Sitter) spacctimc. 
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1 Introduction 

General relativity is no doubt one of the most important and beautiful theories 
of the 20th century. In this theory, Einstein's happiest thought, the equiva- 
lence principle, plays a most essential role. But somewhat ironically it also 
leads us into difficulties in having a proper definition for the energy density of 
gravitating systems, therefore we still do not have a satisfactory description for 
gravitational energy. The modern concept is that gravitational energy should 
be non-local, more precisely quasi-local, i.e., it should be associated with a 
closed two-surface (for a comprehensive review see [I]). Here we consider one 
proposal based on the covariant Hamiltonian formalism fj] wherein the quasi- 
local energy is determined by the Hamiltonian boundary term. For a specific 
spacetimc displacement vector field on the boundary of a region (which can be 
associated with the observer), the quasi- local energy — defined as the value of 
the Hamiltonian boundary term — depends not only on the dynamical values 
of the fields on the boundary but also on the choice of reference values for 
these fields. Thus a principal issue in this formalism is the proper choice of 
reference spacetime for a given observer, which is equivalent to finding a suit- 
able embedding of a neighborhood of the boundary 2-surface into a suitable 
reference space. Of particular interest here is the use of the quasi-local energy 
value to find the embedding variables. 

It is generally accepted that the total energy for an asymptotically flat grav- 
itating system should be non-negative and should vanish only for Minkowski 
space (this is required for stability, see, e.g., [3J; for proofs of this property 
for GR see [HE]). In view of this these properties have also been regarded by 
many as being desirable for a good quasi-local energy [TJEJHHE] ■ 

Following up on one approach used in [TUllllj . recently we proposed an 
encrgy-extremization program [121113] to determine the reference, and found 
that it works well in spherically symmetric spacetimes. We have successfully 
obtained the quasi-local energy for any given displacement vector. However, 
since the program was designed to produce the minimum energy — the program 
chooses one specific embedding as the ground state such that the energy pro- 
duced is the minimum — there is no guarantee of getting non-negative energies. 
We find the stability arguments for positive energy compelling for asymptoti- 
cally flat regions approaching equilibrium. As discussed in some detail in |11[ 
IT2in"3"irPf] , we do not see why a good measure of quasi-local energy for regions 
that are highly dynamic and / or not asymptotically flat must necessarily have 
positive energy. 

Although for our agenda, that is, to find the "best" choice of reference 
for our Hamiltonian boundary term, negative results can be accepted, we still 
wonder under what conditions one can get the energies that satisfy the usual 
criteria, i.e., it should be non- negative and vanishes only for Minkowski (or 
(anti)-de Sitter) spacetime. Our idea is that for any observer in the dynamic 
spacetimc one should pick a comparable reference observer as the ground state. 
Therefore we should minimize the energy under some suitable constraints. 
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There are some common themes and aspects of the work here that are 
similar to those in the recent quasi- local work of Wang and Yau jSUS]. Of 
course their embedding construction is much more general and intricate than 
our spherically symmetric case. In [5] they have used an optimal embedding, a 
strategy that we will use here too. There are also significant differences. First 
and foremost they have a quite general positive energy proof for their energy. 
In the present paper, our spacetime vector is generally allowed to be picked 
by hand, consequently the embedding variables found here usually depend on 
this vector. We have considered here embeddings into a background anti-dc 
Sitter spacetime not just Minkowski. Our embeddings arc not just 2-surface 
isometric embeddings into Minkowski space, they match on the boundary the 
4D dynamic and reference metrics. Also it should be noted that, although 
there are certain similarities, there are also important differences between the 
quasi-local expressions we have developed and those based on the Brown and 
York |15j Hamilton- Jacobi approach, including in particular the expressions 
used by Wang and Yau. 

In this paper we propose a modified energy-extremization program, in 
which we minimize the energy under both the following constraints: 

CI. £ N (^Atf 3 ) = £ N (tf 2 Ad 3 ), 

C2. t?° =0°, N = le = lea, I 2 = -s(N, N), 
where »£n means the Lie derivative along the displacement vector N, and 
means the orthonormal coframe of the dynamic (reference) spacetime. 
The geometric meaning of these constraints is clear. The constraint CI requires 
the time evolution of the two-area for both the dynamic and reference observers 
to be the same on the boundary. It should be kept in mind that throughout 
this paper we always impose the two-surface isometric embedding condition, so 
the two-areas at some constant time are equal too, i.e., i3 2 A rf 3 = r 2 sin 8d8 A 
dej) = d 2 A ?9 3 . Since the quasi-local energy is defined on the two-boundary, 
it is natural to pick a reference observer in this way so we will not compare 
the energy for a expanding two-area with a shrinking one. The constraint 
C2 requires the displacement vector N to be the normal of the spacelike 
hypersurfaces = and d° = 0, and it has the same magnitude for both 
dynamic and reference spacetimes. In other words, we synchronize the two 
clocks and foliate both spacetimes in the same manner such that the flow 
of time is orthogonal to the hypersurface at the boundary 2-surface. These 
constraints determine what we mean by comparable reference observers both 
geometrically and physically. Applied to the spherically symmetric spacetimes, 
the results satisfy the usual criteria. The only exception is that for a small 
region inside the inner horizon of the Reissner-Nordstrom spacetime the energy 
is negative. This negative result is not so strange as this is a rather odd region 
where the gravitational force becomes repulsive. 

In addition to the produced quasi-local energies, it is also very interesting 
and important to investigate the references determined by our energy extrem- 
ization programs. It turns out they are the embeddings that exactly match 
the whole orthonormal frames right on the two-sphere boundary. This result 
motivates us to propose another complementary approach to the whole prob- 
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lem, that is, we propose that the optimal choice of the reference is the one — at 
least for spherically symmetric cases — which matches the 4D metric on the 
boundary. Applying this idea, we find the resultant energy still cannot be 
uniquely determined. This can be understood: in order to get a meaningful 
quasi-local energy we also need to specify a reference observer to compare 
with, so we should further impose some physical conditions. There are two 
obviously physically meaningful choices we can think of. First, it is reason- 
able to require the reference displacement vector to be (proportional to) the 
timclike Killing vector, i.e. N = N T dx- Second, instead of always picking 
the reference timclike Killing vector as our displacement vector, we can also 
pick a comparable reference observer as the ground state. By that we mean 
requiring £n(^ 2 A $ 3 ) = £n{$ 2 A#) on the boundary. As we are going to 
show, these two choices are equivalent to energy extremization programs I 
and II respectively. Therefore the pictures of the quasi-local energy for spher- 
ically symmetric spacetimes are quite complete and satisfactory. On the one 
hand, we know the geometric meaning of the energy extremization programs. 
On the other hand, the most optimized references for spherically symmetric 
spacetimes, i.e. those that match the geometry on the boundary, produce the 
extremum quasi-local energies. 

In the following section we briefly review the quasi-local expression. Then 
apply both energy-extremization programs to spherically symmetric space- 
times in section 3. In section 4 we turn to the isometric matching approach 
and show that the results are compatible with the energy-extremization pro- 
grams. A discussion section concludes this paper. 



2 Quasi-local Expression 

In this section we first briefly sketch the preferred quasi-local energy-momentum 
Hamiltonian boundary term expression for Einstein's gravity theory which was 
derived using the covariant Hamiltonian formalism. The details of the deriva- 
tion and related discussions can be found in [T6l[T7lfT8lll9U20l[2Tj (for some 
additional developments along similar lines see [221123"] ). 

For a theory which is invariant under diffeomorphisms, in particular in- 
finitesimal displacements along some vector field N, the Hamiltonian which 
generates the dynamical evolution of a spatial region along such a vector field 
is given by the integral over the region of a suitable Hamiltonian 3-form which 
takes the form 

U(N) +d»(N). (1) 

As a consequence of local diffcomorphism invariance (i.e., a symmetry for non- 
constant N^), the density part H.^ must be proportional to certain 4-covariant 
field equations, hence it vanishes "on shell". Consequently the value of the 
conserved quantity associated with a local displacement N and a spatial region 
S is determined entirely by a 2-surface integral over the boundary of the 
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region: 

E(N,S):= [ W(N)=* B(N). (2) 

The value is thus quasi-local (depending only on the values of the fields on the 
boundary). For any choice of N this expression defines a conserved quasi- local 
quantity. Several particular boundary terms were identified, each is associated 
with a distinct type of boundary condition. In [2D] a "preferred boundary 
term" (it has a certain covariant property, directly gives the Bondi energy 
flux, and has a positive total energy proof) for GR was identified: 

B(N) = -±-(ALj a p Ai N i 1 J + DpN a AnJ), (3) 

107T 

where = ojI"^ is the connection one-form, r\ a ^ := *($" A i?* 5 ) is a 2- 
form depending on the coframe $ a , A indicates the difference between the 
dynamic and reference values, and Dp is the covariant derivative using the 
reference connection. The reference values can be determined by pullback from 
an embedding of a neighborhood of the boundary into a suitable reference 
space. Now we can use this expression to calculate the gravitational energy in 
general relativity. 



3 Energy-Extremization Program 

In this section, we first briefly review the results (and extend them using 
a generalization from the Schwarzschild to the Reissner-Nordstr6m-(anti)-de 
Sitter (RN-(A)dS) spacetime) obtained in our earlier papers [T2"lll3| . where we 
extremized the energy with respect to the embedding variables without any 
further constraint other than the two-surface isometric embedding condition. 
Then we present our modified program where we extremize the energy under 
some additional suitable constraints. For convenience we call them programs 
I (without constraint) and II (with constraints). Note again that throughout 
this paper we always impose the two-surface isometric embedding condition. 

Program I: 

Extremize the energy with respect to the embedding variables without any 
further constraint. 

We first consider the static spherically symmetric case, more specifically, 
we consider the RN-(A)dS spacetime, 

ds 2 = -Adt 2 + A- X dr 2 + r 2 dn 2 , (4) 

where 

A = 1 - — + % - -r 2 , dn 2 = d9 2 + sin 2 Bd<\> 2 . (5) 
r r 2 3 
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Choose the pure (anti)-dc Sitter spacctime as the reference, 

ds 2 = -BdT 2 + B~ 1 dR 2 + R 2 d(2 2 , (6) 

where B = 1 — ^ R 2 . A legitimate approach for the spherically symmetric case is 
to assume T = T(t, r), R = R(t, r),G = 6,$ = <p along with R := R{t ,r ) 
■ this symmetrically embeds a neighborhood of the two-sphere boundary 
S at (to,ro) into the (A)dS reference such that the two-sphere boundary is 
embedded isometrically. 

Because the spacetime is spherically symmetric, we expect that the quasi- 
local energy can be completely determined by the first term of the expres- 
sion ©. In earlier works, via both explicit calculation and theoretical anal- 
ysis [171124.25,26, 27]. it was found that the second term in our quasi- local 
expression makes important contributions to the value of the angular momen- 
tum and the center-of-mass moment. This term is not expected to make a 
contribution to the energy of spherically symmetric systems, and we have not 
noticed any case where it does. Here in this section we will assume for simplic- 
ity that this holds and will proceed to determine the reference using just the 
first term in our expression. It will turn out that the reference determined by 
this procedure indeed does guarantee that the second term vanishes, so this 
procedure is self-consistent. In the next section we give a further confirmation: 
we will obtain the same results from a different assumption which definitely 
kills the second term. 

Now, using just the first term, from [3]) the quasi-local energy of such a 
spacetime in terms of the displacement vector 

N = N% + N r d r = N T d T + N R d R , (7) 

and embedding variables (T t ,T r , R t , R r ) can be written as 

E = r - (BXN T + AN f R r + A~ 1 N r R t - 2AN t ) , (8) 

where A -1 = T t R r — T r R t . By cxtrcmizing the energy © with respect to the 
four embedding variables we get 

BX 2 T r N R = 0, (9) 
BX 2 T t N R = 0, (10) 
BX 2 T r N T + A _1 A r = 0, (11) 
BX 2 T t N T - AN 1 = 0. (12) 

Note that the first two equations arc equivalent (since we do not want both 
T t and T r to vanish), so we only have three independent equations. Exploiting 



dE 

dE 
df r 
dE 
dR t 
dE 



= 



= 
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these equations we can get 



(13) 

so that the energy obtained is 



Ei = r ( \/-g(N, N)B — AN 1 ) . (14) 



For a static observer at a distance r (note: such observers exist only for 
values of r such that A is positive) 



N static = (15) 
VA 



the energy measured isJZl 



£i(N static ) = 2m Q2 (l (16) 
VB + VA 



For Schwarzschild this reduces to a standard result, I? = r(l — yl — 2m/r), 
first found by Brown and York }15) . 

For a radial geodesic observer who falls initially with velocity vq from a 
constant distance r = a > 2m in Schwarzschild spacetime, i.e., Q = 0,A = 0, 



the energy measured according to our program is 



E ^o)=r(l-^^ T ^y, (18) 

(values which were also found by other techniques, see [2"Bll2"5] ). When the 
initial velocity vo is less, equal, or greater than the escape velocity, y / 2m/a, 
the energy is positive, zero, or negative, respectively. One interesting fact is 
that the sign of the scalar curvature of the hypersurface orthogonal to this 
displacement vector, N geo , is the same as the energy. This feature is very 
much like the cosmology case which we are going to see later. 

It is obvious that the energy-extremization equations (|5r ll2[) cannot uniquely 
determine the reference, since only three of them are independent. However, 
by further imposing the condition g(N, N) = ^(N,N) on the boundary the 
reference determined is unique and kills the second term of the expression ((3]) . 



This is meaningful only for ranges of r where B and A are positive. A similar remark 
applies to many of our later expressions. We will often not explicitly mention the appropriate 
domain of validity; it can easily be determined by a detailed examination of the expressions. 
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Program II: 

Extrcmizc the energy with respect to the embedding variables under both the 
following constraints. 

CI. £ N (ti 2 _ A?? 3 ) = £ N (^ 2 Ai? 3 ), 

C2. #° = N = le = le , I 2 = -g(N, N), 
where £n means the Lie derivative along the displacement vector N. 

From constraint CI we get 

£ N (t? 2 A i? 3 ) = 2r sin 6W r rf6» A d(f>, 
£ N (tf 2 A i9 3 ) = 2rsm6N R d6 A d</>, 

i?N(t? 2 Atf 3 ) = £ N (tf 2 Atf 3 ) N r = N R , (19) 

where the useful formula »£n = zn^ + g^n is used. Note that these constraints 
are imposed only on the boundary, where (because of isometric matching) 
R = r. Then from (fT9|) we get the relation between R t and R r , 

N t R t + N r R r = N R = N r ^ = _(i_ jRr ). (20) 

In order to apply constraint C2, we first write the coframes in the form 

if = a t dt + a r dr, i? 1 = b t dt + b r dr, 
§° = a T dT + a R dR, tf 1 = b T dT + b R dR. (21) 

This is quite general since we only consider spherically symmetric cases here; 
what is being allowed are radial boosts of the obvious coframes associated 
with ((4J and (j6|) . Note that one can also express the displacement vector using 
the orthonormal frames 

N = N t d t + N r d r = N T d T + N R d R 

= N°e + N 1 e 1 = N°e + N l e x . (22) 

Using ((21]) and ([22j we get: 

N° = #°(N) = a t N l + a r N r , 
N 1 = ^(N) = b t N l +b r N r , 
N° = d°(N) = a T N T + a R N R , 
N 1 = ^(N) = b T N T + b R N R . 

Together with 

-Adt 2 + A- l dr 2 = -(i? ) 2 + (d 1 ) 2 
=> -a 2 t +b 2 t =-A, -a 2 + b 2 r = A- 1 , -a t a r + b t b r = 0, (23) 

and 

-BdT 2 + B- l dR 2 = -(i9 ) 2 + (d 1 ) 2 
=S> -ciT + b 2 T = -B, —a% + b 2 R = B' 1 , -a T a R + b T b R = 0, (24) 
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one can work out the following relations: 

_ AN ^ t N x N r -A~ 1 N°N r ± N 1 ^ 

at - -,g(N,N) ' -,g(N,N) ; 

AN 1 ^ t N°N r -A~ 1 N 1 N r ± N ^ 

-ff(N,N) ' r ~ -9(N,N) ' 



_ BN°N T T N J N R —B~ 1 N°N R ± iV 1 ^ 

_bn 1 n t tN°n r -b- 1 n i n r ±n°n t 

T ~ -fl(N,N) ' R ~ -g(N,N) ' 



(25) 



In order to preserve the orientation, one should choose the upper sign. Now 
by applying constraint C2 we get 

i9° = a t dt + a r dr = a T dT + a R dR = t?° 
=> at = a T T t + aRRt, a r = axT r + anR r , (26) 

and 

N = le = le => iV 1 = iV 1 = 0. (27) 
By applying (|25M^T|) we get the following relations 

T t = ^ ¥ (AN t + B- 1 N R R t ), 

T r = -^{-A- 1 ^ + B- 1 N R R r ). (28) 

Together with (|20p then there is now only one independent variable, which we 
choose to be R r . Extrcmizing the energy with respect to R r gives 

Rr =^N ) ( BNTNt ~ A ~ 1 ^ 2 y 

Note that by straightforward calculation one can see that the determined 
reference kills the second term of The energy obtained is 

En = r (^y/-g{N,N)B + (N r ) 2 - AN*} 
-g(N,N)(2m-Q 2 /r) 



y/-g(N, N)B + (N r ) 2 + AN* 



(30) 



One can see that for any timelike displacement vector this energy is always 
non-negative except for a small region inside the inner horizon of the RN 
black hole, r < Q 2 /2m. The transition is exactly at the turning point where 
the gravitational force becomes repulsive. For the case m = Q = the energy 
vanishes for all displacement vectors, just as it should, whereas in the case 
m = Q = 0, (TH| vanishes only when N r = 0. 

For a static observer, the result is the same as program I, 

£n(N static ) = £i(N static ). (31) 
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For the radial geodesic observer in Schwarzschild spacctime the energy mea- 
sured is 

9m 

^i(N g co) = - j==, (32) 

which is positive for any initial velocity ^o- 

For the dynamic case, these programs work well too for the Friedman- 
Lemaitre-Robertson- Walker (FLRW) metric. The FLRW metric in the stan- 
dard spherical coordinate is 

ds 2 = -dt 2 + A 2 dr 2 + a 2 {t)df} 2 , A= ^ . (33) 

v 1 — kr 2 

Apply program I, the resultant energy obtained from expression ([3]) is 



Ei = ar ( V-#(N, N) - i _1 aiV* - AarN r ) , (34) 



where a = da/ dt. 

For a comoving observer this gives a value that was found earlier in |30j . 

Ej(d t ) = ^ (35) 

V 1 1 + v/1 - kr 2 V ' 

which is positive, zero, and negative for k = +1,0,-1 respectively. An in- 
teresting geometrically preferred observer is the one preserving the area. The 
associated evolution vector is the dual mean curvature vector, known as the 
Kodama vector. In this case we have 



y/1 - kr 2 ar \fi - kr 2 

JN Kodama - , -99 °^ /-, , 9 =^r^° T - \ 6K) ) 

yl - kr z — a z r z a yl- kr z — a z r z 
The corresponding energy is 

p, N ^_ ar 3 (fc + d 2 ) 

El(S\ Kodama ) - — , , ==, (37 ) 

1 + vl — kt — a z r z 
which is positive when the Friedman equation, k + a 2 = ^ pa 2 , is satisfied: 

F rN ¥KH 3 _ 2m(f) 

(,JM Kodama j — , — , = , (OO) 



3 



. p ( ar )2 _ ?m(f) 



where m(f) = ^-pr 3 and f = ar. 



On the other hand, applying program II we get 

^Il(lNKodama) - ^ - ( . _ - + ^^3^^ + ^-S^r ' ^ 



where Z = — g(N, N). This energy is always positive for any timelike displacc- 

3 



ment vector N when the Friedman equation, k + a 2 = %-pa 2 , is satisfied. Note 
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that it vanishes for the Milne universe, i.e. k = — 1, a = t, which is diffcomor- 
phic to Minkowski spacetimc. 

For a comoving observer, the energy measured is 

E , 9)= ar 3 (k + a 2 ) feV ( ) 

{ t} y/1 + a 2 r- 2 + VI - kr 2 y/l + a 2 r 2 + y/l- kr 2 ' V ; 

For a Kodama observer the result, after substituting TV*, N r from ([55)1 . works 
out to be the same as was found using program I: 

-Ell(NKodama) = E\ (Nxodama) • (41) 



4 An Alternative Approach 

We have seen that the proposed optimal programs give reasonable quasi-local 
energies: it is natural then to ask what are the determined references. For 
simplicity let us consider the Schwarzschild geometry in the standard spherical 
coordinate and the unit future timelikc displacement vector N. The embedding 
variables determined by program I are |121ll3j 

T t = AN*, T r = -A _1 iV r , R t = -N r , R r = N*. (42) 

It is not very difficult to see that for N to be the unit normal of the constant 
time hypersurface, i.e., 

N = e = —j=dt = d T = e 
VA 

=> JV* = -=, N r = 

T t = VA, T r = 0, Rt = 0, R r = -^=, 

VA 

the embedding isometrically matches the complete 4D orthonormal frame at 
the two-sphere boundary, this means the additional restrictions 

=dT = T t dt + T r dr = VJdt = 0°, 

d 1 =dR = R t dt + R r dr = -^=dr = d 1 . 

VA 

This is not extremely surprising since we are only considering spherically 
symmetric spacetimes here. So we propose an alternative approach here: the 
optimal choice of the reference should be the one which isometrically matches 
the 4D orthonormal frame at the two-sphere boundary, i.e., on the boundary 
we require 

i9° = 0° ^> a t dt + a r dr = a T dT + a R dR, 

i? 1 = 1 =S> b t dt + b r dr = b T dT + b R dR, (43) 
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which implies 



and 



eo = e , ei=ci, N° = iV 5 , 

iV 1 = N 1 , 5 (N,N) = 5(N,N), (44) 



Tt = a t &R - a H 6 t , i? t = -a t b T + a T b t , 
T r = a r b R — a R b r , R r = —a r bT + arb r . (45) 

For these embeddings the second term of vanishes. By ([25"]) we get the 
following results: 

X 1 = T t R r — T r R t 

= {a t bR - a R b t ){~a r b t + a T b r ) - {a r b R - a R b r )(-a t b T + a T b t ) 
= 1, 

and 

AN l R r + A~ x N r R t = -(AN^r + A" 1 iV r a t )6 T + (AN% + A _1 iV r 6 t )a T 
= -N x b T + N°a T = BN T , 

so the energy obtained, for the RN-(A)dS metric, is 

E iso = r{BN T - AN') 



,Jl 2 B + (N R ) 2 - ^l 2 A+(N r )^ 

l 2 (2m - Q 2 /r) + r{(N R ) 2 - (N r ) 2 ) 
^l 2 B + (N^+^PA+(N r ) 2 



(46) 



where I 2 = -g(N,N). 

For any given displacement vector in the dynamic spacetime there is still 
one unknown in this energy expression, N T (or N R ). Physically to determine 
it means to pick a corresponding observer in the reference spacetime as the 
ground state. There are two obvious and physically meaningful choices we can 
think of, and, as we are going to show, they are equivalent to the energy- 
cxtrcmization programs I and II respectively. First, it is reasonable to require 
the reference displacement vector to be (proportional to) the timelike Killing 
vector, i.e., N = N T dT- Then (|4"6"]l reduces to 



E iso i (N) = r v/-g(N,N)B -AN'), (47) 



which is equivalent to (fT4")l . Second, it is also reasonable to pick a comparable 
reference observer as the ground state. By that we mean requiring, on the 
boundary, 

4(^Ai? 3 )=4(^AiJ 3 ), (48) 
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where £jsi means the Lie derivative along the displacement vector N. This 
condition implies 



N R = N r => BN T = y/—g(N, N)B + (N r ) 2 (49) 
so that (l46l) reduces to 



£isoii (N) = r ( v / -.9(N,N)B+(iV-) 2 - AN*) , (50) 



which is equivalent to ((30)) . 

So one can see that this isometric matching approach gives results equiva- 
lent to the energy extremization programs. These two methods can be seen as 
complementary to each other. It helps us understand the geometric meaning 
of the energy-extremization program, at least for the spherically symmetric 
spacetimes, on the one hand. On the other hand, when we choose the refer- 
ence that isometrically matches the 4D metric, the resultant energy would be 
the extremum. 

Note that the quantities in the dynamic spacetime are known, and N° = 
N° and N 1 = N 1 after the exact matching (|g5]t. So from (gSJ the coefficients 
of the reference orthonormal coframes (o-t, i.r, &T) & r) are functions of the 
embedding variables (T t ,T r , R t , R r ), since 

N T = N t Tt + N r T ^ N R = N t Rt + N r R ^ ^ 

Then one can read (|24p as three constraints on (T t ,T r , R t , R r ), so that now 
there is only one independent variable in (|46p . It is natural then to expect that 
one can reproduce E\ if one extremizes the energy (|4*5| with respect to the sin- 
gle variable. Technically instead of the embedding variables (T t , T r , R t , R r ), it 
is more convenient to vary (|46)) with respect to the coefficients (ar, an, br, &r)- 
For simplicity we normalize the displacement vector, 

ff (N,N) = -l = 3(N,N). (52) 

From the equations for or and 6t in (|25p we get 

NT = N^NHr 
B 

Now by extremizing (|46[) with respect to ax we get 

= N o_ N i^L (54 ) 
oa T 

which together with (|24|) gives 
dbx o-t 

N 1 N° 
a T = y/BN°, b T = VBN 1 , a R = -=, b R = — , 

v B v B 



N T = -L. (55) 
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Table 1 The energies for the Reissner-Nordstr6m-(anti)-de Sitter and Friedman-Lcmaitrc- 
Robertson- Walker spacetimes by different approaches. 



energy for RN-(A)dS energy for FLRW 



iso 


r(BN T - AN 1 ) 






ar(N T - 2^ - 


Adr-ZV") 






isol 


r(V-s(N,N)B- 


A7V 4 ) 




ar (v^2 _ a^i _ 


AarN r ) 






isoll 


r(V- 9 (N,N)B + 


(A?r)2 


— AN*) 




f aW) 2 - 


aJV* 
A 


- AarN r ) 


I 


r(V-<?(N,N)B- 


AN*) 




ar (v^2 - - 


AdrAT r ) 






II 


r(V~9(N,N)B + 


(N r ) 2 


-AN*) 


ar( % /i 2 + (drAT t - 


f aiV) 2 - 


aiV* 
A 


- AarN r ) 




A=1 _2rn + & 

r r z 

B = 1- |r 2 


3 ' 




A - " 

Vl-fer 2 

Z 2 = — g(N, N) 









iso means given (7V i ,7V r ) then matching the orthonormal frames, 
isol means iso with the restriction N R = 0. 
isoll means iso with the restriction CI. 

I means given (N t ,N r ), extremize the energy with no constraint. 

II means given (N t ,N r ), extremize the energy under the constraints CI and C2. 



So the energy obtained is 

£isocx = r(\/S - AN') = E 1 = E isol (56) 

as expected. 

Applying this isometric matching approach to the FLRW spacetime we get 



E iso = ar (N- - yfr=k*N* - - ^ j . (57) 

If wc require N = N T dx then (fSTf is equal to On the other hand, by 

requiring pg|) . which implies N R = arN 1 + aN r , the energy is the same as 
program II. If we extremize (|57[) with respect to the only one independent 
variable, the result is the same as that determined by program I (|34j) . 

Since we have explored so many different approaches to determine the 
reference and quasi-local energy, it may be suitable here to summarize these 
various approaches to quasi-local energies for spherically symmetric spacetimes 
by the Table [TJ 



5 Discussion 



As mentioned in the introduction, a principal issue in the covariant Hamil- 
tonian formalism is the proper choice of reference spacetime for a given ob- 
server. When we embed the two-sphere boundary in the dynamic spacetime 
into the reference, the question becomes how to determine the embedding. In 
our earlier papers [T2UT5] we imposed the isometric embedding condition and 
minimized the energy with respect to the embedding variables with no further 
constraint, so the energy produced is the absolute minimum of all possible 
embeddings for any given displacement vector. For all different observers in 
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the dynamic spacetime the program has chosen the same reference observer 
to be the ground state, so for some observers with extreme motion the energy 
measured can be extreme. In this paper we have proposed a modified program 
in which we minimize the energy under some additional suitable constraints. 
The main idea is to force the program to pick a comparable reference observer 
while minimizing the energy. Then no matter how extreme the motion is of 
an observer in the dynamic spacetime, he will be taking a comparable refer- 
ence observer as the ground state. Then the results satisfy the usual criteria. 
The only exception is that for a small region inside the inner horizon of the 
Rcissner-Nordstrom spacetime the energy is negative. However, it is rather an 
odd region where the gravitational force becomes repulsive. 

We have also shown that the 4D isometric matching approach produces 
compatible results. These two methods, energy-extremization and 4D isomet- 
ric matching, can be seen as complementary to each other. It helps us under- 
stand the geometric meaning of the energy-extremization program, at least 
for the spherically symmetric spacetimes, on one hand. On the other hand, 
the most optimized references for spherically symmetric spacetimes, i.e., those 
that match the geometry on the boundary, produce the extremum quasi-local 
energies. We believe that this spherical case is the main test case; it shows that 
our programs have promise as universal approaches for determining the refer- 
ence needed for the covariant Hamiltonian boundary term quasi-local energy 
for general spacetimes. 
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